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Abstract 



To generate small neutrino masses radiatively, the Zee-model introduces two 
Higgs doublets and one weak-singlet charged Higgs boson to its Higgs sector. 
From analyzing the renormalization group equations, we determine the possi- 
ble range of the lightest CP-even Higgs boson (h) mass and the Higgs boson 
self-couplings as a function of the cut-off scale beyond which either some of 
the coupling constants are strong enough to invalidate the perturbative anal- 
ysis or the stability of the electroweak vacuum is no longer guaranteed. Using 
the results obtained from the above analysis, we find that the singlet charged 
Higgs boson can significantly modify the partial decay width of /i ^ 77 via 
radiative corrections, and its collider phenomenology can also be drastically 
different from that of the charged Higgs bosons in the usual two-Higgs-doublet 
models. 

PACS number(s): 13.15. +g, 14.80.Cp. 
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I. INTRODUCTION 



Prom the atmospheric and solar neutrino data, there is increasing evidence for neutrino 
oscillations [1]. If this is a correct interpretation, the Standard Model (SM) has to be 
extended to incorporate the small masses of the neutrinos suggested by data. There have 
been several ideas proposed in literature to generate small neutrino masses. The Zee-model 
is one of such attempts [2-6]. In this model, all flavor neutrinos are massless at the tree 
level, and their small masses are induced radiatively through one-loop diagrams. For such a 
mass-generation mechanism to work, it is necessary to extend the Higgs sector of the SM to 
contain at least two weak-doublet fields and one weak-singlet charged scalar field. Although 
some studies have been done to examine the interaction of the leptons and the Higgs bosons 
in the Zee- model [7], the scalar (Higgs) sector of the model remains unexplored in detail. 
In this paper we study the Higgs sector of the Zee-model to clarify its impact on the Higgs 
search experiments, either at the CERN LEP-II, the Run-II of the Fermilab Tevatron, the 
CERN Large Hadron CoUider (LHC), or future hnear coUiders (LC's). 

Experimental search for the Higgs boson has been continued at the CERN LEP and the 
Fermilab Tevatron experiments. In the LEP-II experiments, the Higgs boson with the mass 
less than about 110 GeV has been excluded if its production cross section and decay modes 
are similar to those of the SM Higgs boson [8]. Run-II of the Tevatron can be sensitive to 
a SM-like Higgs boson with the mass up to about 180 GeV, provided that the integrated 
luminosity of the collider is large enough (about 30fb~^) [9]. Furthermore, the primary goal 
of the CERN LHC experiments is to guarantee the discovery of a SM-hke Higgs boson for its 
mass as large as about 1 TeV [10], which is the upper bound of the SM Higgs boson mass. 
(For a Higgs boson mass beyond this value, the SM is no longer a consistent low energy 
theory.) 

When the Higgs boson is discovered, its mass and various decay properties will be mea- 
sured to test the SM and to distinguish models of new physics at high energy scales. For 
example, the allowed mass range of the hghtest CP-even Higgs boson (h) can be determined 
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by demanding the considered theory to be a vahd effective theory all the way up to some 
cut-off energy scale (A). For A — GeV (i.e., the Planck scale), the lower and upper 
bounds of the SM Higgs boson masses are 137 GeV and 175 GeV, respectively [11]. The 
Higgs mass bounds for the two-Higgs-doublet-model (THDM) were also investigated [12,13] 
with and without including the soft-breaking term with respect to the discrete symmetry 
that protects the natural flavor conservation. It was found in Ref. [13] that the lower bound 
of the hghtest CP-even Higgs boson is about 100 GeV in the decouphng regime where only 
one neutral Higgs boson is light as compared to the other physical states of Higgs bosons. 

The Higgs sector of the Zee-model is similar to that of the THDM except for the existence 
of an additional weak-singlet charged Higgs field, so that the physical scalar bosons include 
two CP-even, one CP-odd and two pairs of charged Higgs bosons. In this paper, we shall 
first determine the upper and lower bounds for the lightest CP-even Higgs boson mass (m/j) 
as a function of the cut-off scale A of the Zee-model, using renormalization group equations 
(RGB's). ^ We show that the upper and lower mass bounds for h are almost the same as those 
in the THDM. We also study the possible range of the Higgs-boson self-coupling constants 
at the electroweak scale as a function of A. By using these results, we examine effects of the 
additional loop contribution of the singlet charged Higgs boson to the partial decay width of 
h 77. We show that, by taking A = 10^^ GeV, the deviation of the decay width from the 
SM prediction can be about —20% or nearly -|-10% for rrih between 125 GeV and 140 GeV 
when the mass of the isospin singlet charged Higgs boson is taken to be around 100 GeV. 
The magnitude of the deviation becomes larger for lower cutoff scales and smaller masses of 
the singlet charged Higgs boson. If we choose A = 10^ GeV and the singlet charged Higgs 
boson mass to be 100 GeV, the positive deviation can be greater than -|-30% (-1-40%) for 
rrih = 125 GeV (140 GeV). Such a deviation from the SM prediction could be tested at the 



For the model with see-saw mechanism for neutrino mass generation the Higgs mass bound has 
been studied as a function of cut-off scale in Ref. [14] . 
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LHC, the e+e~ LC and the 77 option of LC [15-17]. We also discuss phenomenology of the 
singlet charged Higgs boson at present and future collider experiments, which is found to be 
completely different from that of the ordinary THDM-like charged Higgs bosons. To detect 
such a charged Higgs boson at LEP-II experiments, experimentalists have to search for their 
data sample with or /j,"^ plus missing energy, in contrast to the usual detection channels: 
either ru or cs decay modes. 

This paper is organized as follows. In Sec. II, we introduce the Higgs sector of the Zee- 
model and review the neutrino masses and mixings in this model which are consistent with 
the atmospheric and solar neutrino observations. Numerical results on the possible range 
of the mass and coupling constants of the Higgs bosons are given in Sec. HI. In Sec IV, we 
discuss the one-loop effect of the extra-Higgs bosons in the Zee-model to the partial decay 
width of /i — > 77 and its impacts on the neutral Higgs-boson search at high-energy colliders. 
The phenomenology of the charged Higgs boson that comes from the additional singlet field 
is discussed in Sec. V. In Sec. VI, we present additional discussions and conclusion. Relevant 
RGE's for the Zee-model are given in the Appendix. 



II. ZEE-MODEL 

To generate small neutrino mass radiatively, the Zee-model contains a SU{2)l singlet 
charged scalar field a;~, in addition to two SU{2)l doublet fields 0i, and 02- The Zee-model 
Lagrangian is written as: 

— ^kin + ^lliM + ^Yukawa — V{(f>l, 02, <^ ) , (l) 

where 

2 

+U;rD,l,+ie^rD,e,+ ^ ^Pf,, (2) 

o=5C/(3),5!7(2),!7(l) ^ 



where i,j {— 1, 2, 3) are the generation indices, and 



y(0i,02,t^ ) ^mj\(f)if + ml\(f)2f + ml 



+ |0i|^ + |02|^ + A3 \(f)2f 

+A4 

+ai CO \0i\' + (J2 CO \02\' + -o-^ OJ 



2 

2 A5 
2 



|0l| +0-20; |02| + ^0"3 



(4) 



In the above equations, ql is the left-handed quark doublet with an implicit generation index 
while ur and denote the right-handed singlet quarks. Similarly, II and Cr denote the left- 
handed and right-handed leptons in three generations. The charge conjugation of a fermion 



field is defined as ■0*^ = C^l) , where C is the charge conjugation matrix {C~^^^C = —7'*^) 
with the super index T indicating the transpose of a matrix. Also, (f)^ = 



and 



y Vm J 



(f)jn = (ir2) (j)^ with m = 1,2. Without loss of generality, we have taken the anti-symmetric 
matrix f^j and the coupling /i to be real in the equations (3) and (4). In order to suppress 
flavor changing neutral current (FCNC) at the tree level, a discrete symmetry, with 0i (f>i, 
4>2 —4>2, cu~^ —>■ +co~^, is imposed to the Higgs sector of the Lagrangian, which is only 
broken softly by the m| term and the /i term. Under the discrete symmetry there are two 
possible Yukawa-interactions; that is, for type-I 

'-'Yukawa— I 

and for type-II, 



■'Yukawa— II 



dR. {voycKM).. 'PiQl^ + {yu)ii 4>2(1l, + {yE)ii 'Pik, + h.c. , (6) 



where y^,y^,yj^ are diagonal Yukawa matrices and Vckm is the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. Later, we shall only keep the top Yukawa coupling constants yt = {yu)^^ 
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in our numerical evaluation of the RGE's^. In that case, there is no difference between the 
Yukawa couplings of the type-I and the type-II models. Finally, for simplicity, we assume 
that all Xi and are real parameters. 

Let us now discuss the Higgs sector. The SU{2)lxU{1)y symmetry is broken to C/(l)em 
by and (02)) the vacuum expectation values of 0i and (p2- (They are assumed to be 
real so that there is no spontaneous CP violation.) The number of physical Higgs bosons 
are two CP-even Higgs bosons {H,h), one CP-odd Higgs boson (^4) and two pairs of charged 
Higgs boson {Sf, Sf). We take a convention of mn > ruh and > m^j. In the basis 

/ ,0\ 

where two Higgs doublets are rotated by the angle /3, with tan/3 = mass matrices 



for the physical states of Higgs bosons are given by 



M 



N- 



( Ai cos"^ P+\2 sin^ /3 + 1 sin^ 2/?) ( A2 sin^ /? - Ai cos^ /? + A cos 2/?) 
(A2 sin2 /3 - Ai cos2 /3 + A cos 2(3) '^v^ + (Ai + A2 - 2A) 2^M„2 



for CP-even Higgs bosons. 



(7) 



Ml = M^- Asv^ 



(8) 



for CP-odd Higgs boson, and 



2 ^ 



V2 



'V2 



rrin 



(9) 



^0+ (f cos2/3+ f sin2/3) v' 
for charged Higgs bosons. Here, A = A3 -|- A4 -|- A5 and = m|/ sin /3 cos /3. The vacuum 



expectation value v (~ 246 GeV) is equal to \/2^ + (02)^- Mass eigenstates for the 
CP-even and the charged Higgs bosons are obtained by diagonahzing the mass matrices (7) 
and (9), respectively. The original Higgs boson fields, 0i, 02, i can be expressed in terms 
of the physical states and the Nambu-Goldstone modes and G^) as 



^Our analyses will thus be valid in the cases where the effect of the bottom Yukawa coupling is 
sufficiently small; i.e. in the region of not too large tan (3. 
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4>l = ^= (v cos (3 + H cos a — hsina + i{G° cos f3 — Asm 
(f)i — G~ cos (3 — {Si cos x~ ^2 sin x) sin 

02 = sin /9 + if sin a + /i cos a + i{G^ sin (3 + A cos 

(f)2 — G~ sin P + {Si cos x ~ ^2 sin x) cos /3, 
a;- = 5-- sinx + -S^ cosx , 



(10) 

(11) 
(12) 

(13) 

(14) 



where the angle a and x ^-re defined from the matrices which diagonahze the 2x2 matrices 
and M|, respectively. Namely, we have 



cos(a — /5) sin(a — /3) 



^ — sin(Q; — (3) cos{a — (3) j 



N 



cos{a — f3) — sin(a — /5) 



sin(Q; — (3) cos{a — 13) j 



m 



H 







cos X ~ sm X 



\ 



cos X sm X ^ ^2 
y-sinxcosxy \^smx cosx 

where m^j > m\ and > . The mixing angles a and x then satisfy 

M2 - (As + A4 + A5) 



ml, 



Si 

"^&y 



(15) 



(16) 



tan 2a — 



tan2x 



M2-(AiCOs2/?-A2sin2/?)^ 
—y/2fj,v 



tan 2/3, 



(17) 
(18) 



M2 -mg - (A4 + A5 + (7iCos2/3 + (72sin2/3) f ' 

which show that a and x approaches to /3 — and zero, respectively^, when is much 
greater than v^, /x^ and mg; i.e., in the decoupling regime. In this limit, the massive Higgs 
bosons from the extra weak-doublet are very heavy due to the large M so that they are 
decoupled from the low energy observable. 

Although neutrinos in this model are massless at the tree level, the loop diagrams in- 
volving charged Higgs bosons, as shown in Fig. 1, can generate the Majorana mass terms 
for all three-flavors of neutrinos. It was shown [2] that at the one-loop order, the neutrino 
mass matrix, defined in the basis where the charged lepton Yukawa-coupling constants are 



Recall that we assumed 171^ > my^. 
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u> 



3 



FIG. 1. A representative diagram that generates the neutrino mass. For type-I, i = 1, j = 2, 
and for type-II, i = 2, j = 1. 

diagonal in the lepton flavor space, is real and symmetric with vanishing diagonal elements. 
More explicitly, we have 

mi2 mi3 



"7.12 777.23 

^mia m23 



(19) 



with 



fij{ml. - ml) II cot (3 



167r2 



■In"^^^ 



- "^52 "'52 

where (i = 1,2,3) is the charged lepton mass for type-I. For type-II, cot /? should be 
replaced by tan /3. Note that Eq. (20) is vahd for rug.^ rriey 

The phenomenological analysis of the above mass matrix was given in Ref. [4,5]. It was 
concluded that, in the Zee- model, the bi-maximal mixing solution is the only possibility to 
reconcile the atmospheric and the solar neutrino data. Here we give a brief summary of their 
results, for completeness. Let us denote the three eigenvalues for the neutrino mass matrix, 
cf. Eq. (19), and 7771,3, which satisfy 777,^^ -|- 777,^2 -|- 7711,3 — 0- The possible pattern 

of the neutrino mass spectrum which is allowed in the Zee-model is |77i,,J ~ liTi^^l ^ |777j,3|, 

0(10-3) 

eV^ from the atmospheric neutrino data, and Aml^i^^ — 0(10"^) eV^ (MSW large angle 
solution) or O{10~^°) eV^ (vacuum oscillation solution) from the solar neutrino data^. Thus, 



777 



2 ' 



(20) 



with ml^-ml.^ ~ 777 "''^Is = ^^^Itm^ ^^^^ 



^1, - 



^"^Lar> where A7r72^^ 



^Due to the structure of the mass matrix, cf. Eq. (19), only the hierarchy pattern \mi,^ | ~ IrUiy^ | 
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we have \m^^ \ ~ \mi,r^\ ~ ^J^rr^^ {rrii,^ ~ —fn^j^) and |m; 
form of the neutrino mass matrix is given by 

/ 



V3 I 



Am2 , 

solar 











2 


2 



The approximate 







(21) 



where the upper (lower) sign corresponds to rriui < (> 0) case, and the corresponding 
Maki-Nakagawa-Sakata (MNS) matrix [18] which diagonalizes the neutrino mass matrix is 



U = 



+ 



'V2\ 







+ 











^ \/2y I mi/j I + 1 I 



1 

V21 
1 

V2I 











+ 






/ 


m„2 





'"2 







1 

V2 

1 

V2 



(22) 



In the above equations, we took the hmiting case where t/13 = and C/32 = C/23 = ^■ 
From Eqs. (20) and (21), we obtain 



12 



/l3 

/l3 



^ ~ 3 X 10^, 
mi 



f'. 



23 



V2A 



(23) 
(24) 



10^ (for the MSW large angle solution) 
10^ (for the vacuum oscillation solution) 

Therefore, the magnitudes of the three couphng constants should satisfy the relation I/12I S> 
I/13I ^ 1/23 1- This hierarchy among the couplings fij is crucial for our later discussion on 
the phenomenology of the singlet charged Higgs bosons. 

For a given value of the parameters ms-^-, Tns2, tan/? and fi, the coupling constants fij can 
be calculated from Eq. (20). For example, for ms^ — 500 GeV, ms2 — 100 GeV, tan/3 = 1, 
fjL — 100 GeV and mi2 = 3 x 10~^eV, we obtain I/12I ~ 3 x 10~^. As in this example, when 



Imj^gl, rather than \mi,-^ \ ~ |m^.2l ^ l^i^sl) is realized in the Zee-model [4,5] 
^This limit corresponds to 62 = j and ^3 = in the notation of Ref. [18]. 
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Si is rather heavy and the hghter charged Higgs boson S2 is almost a weak singlet, i.e. the 
mixing angle x approaches to zero, it is unlikely that there are observable effects to the low 
energy data [7]; e.g., the muon life-time, the universality of tau decay into electron or muon, 
the rare decay oi /i ^ 67, the universality of VF-boson decay into electron, muon or tau, 
and the decay width of Z boson. When \fij\ are small, we do not expect a large rate in the 
lepton flavor violation decay of a light neutral Higgs boson, such as /i — > /I'^e^ (the largest 
one), h — > e^r^, or /i — > h'^t^ (the smallest one). On the contrary, as to be discussed in 
Section IV, the decay width of ^ 77 can significantly deviate from the SM value. 

Finally, the phenomenological constraints on /12 were derived in Ref. [6]. From the 
consistency of the muon decay rate and electroweak precision test it was found that 

43<7x IQ-^Gf , (25) 

M 

where Gp is the Fermi constant, and 

1 sin^Y cos^Y 
=2 = + 26 

This means that the fij cannot be 0(1) unless the charged Higgs boson masses are at the 
order of 10 TeV. 



III. HIGGS BOSON MASS AND COUPLINGS THROUGH RGE'S 

In this section, we determine the bounds on the mass of the lightest CP-even Higgs boson 
as a function of the cut-off scale of the Zee-model by analyzing the set of renormalization 
group equations (RGE's). We also study the allowed ranges of the coupling constants, 
especially ui and 02 in Eq. (4). In Sec. IV, they will be used to evaluate how much the 
partial decay width of /i — > 77 can deviate from its SM value due to the one- loop contribution 
from the singlet charged Higgs boson. 

The mass bounds are determined in the following manner. For each set of parameters 
defined at the electroweak scale, the running coupling constants are calculated numerically 



11 



through RGE's at the one-loop level. We require that all the dimensionless couphng con- 
stants do not blow up below a given cut-off scale A, and the coupling constants satisfy the 
vacuum stability condition. We vary the input parameters at the electroweak scale and 
determine the possible range of the hghtest CP-even Higgs boson mass as a function of A. 
In a similar manner, we also study the allowed ranges of various Higgs boson self-coupling 
constants at the electroweak scale as well as a function of the lightest CP-even Higgs boson 
mass. 

We derived the one-loop RGE's for the Zee-model, and hsted them in the Appendix for 
reference. For simplicity, in the RGE's, we have neglected all the Yukawa coupling constants 
iUu, Ud-i He) but the top Yukawa coupling y^.^ Although we kept the new coupling constants 
fij in the RGE's listed in the Appendix, we have neglected fij in the numerical calculation. 
This is because the magnitudes of these coupling constants are numerically too small to 
affect the final results unless the singlet-charged scalar-boson mass is larger than a few TeV 
[ cf . Eq. (25) ] . The dimensionless couphng constants relevant to our numerical analysis are 
the three gauge- coupling constants gi, g2, g3, the top Yukawa- coupling constant yt, and eight 
scalar self-coupling constants, Xi {i — 1 — 5) and Uj {i — 1 — 3). There are five dimensionful 
parameters in the Higgs potential, namely mf, m|, m|, mg and Instead of , m|, m|, 
we take v, tan (3, and = m|/ sin (3 cos as independent parameters, where v (~ 246 
GeV) characterizes the weak scale and M the soft-breaking scale of the discrete symmetry. 

In the actual numerical calculation we first fix tan /3 and M. For a given mass {nih) of 
the lightest CP-even Higgs boson, we solve one of the Aj, which is chosen to be A3 here, 
in terms of other Aj. We then numerically evaluate all dimensionless coupling constants 
according to the RGE's. From ruh to M we use the SM RGE's, which are matched to the 



^In the model with the type-II Yukawa interaction, the bottom-quark Yukawa interaction can 
become important for a large tan p. 
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Zee-model RGE's at the soft-breaking scale M. ^ 

We requires the following two conditions to be satisfied for each scale Q up to a given 
cut-off scale A. 

1. Applicability of the perturbation theory implies 

A,(g)<87r, a,(Q)<87r, y^,{Q) < At: . (27) 

2. The vacuum stability conditions must be satisfied. The requirement that quartic 
coupling terms of the scalar potential do not have a negative coefficient in any direction 
leads to the following conditions at each renormalization scale Q: 



(a) 



Ai(g)>0, A2(g)>0, a3(g)>0. (28) 



(b) 



.,(g) + .MiM-M>o, (29) 



,,(g) + . M%M>o, (30) 



A(g) + ^Ai(g)A2(Q) >0, (31) 
where A(g) = A3(g) + min (0, A4(g) + A5(g), A4(g) - A5(g)). 



The parameter tuq and /i arc only relevant to the charged scalar mass matrix. In principle, our 
numerical results also depend on these parameters through the renormalization of various coupling 
constants from the scale of ruh to the charged scalar mass. Since these effects are expected to be 
small, we calculate the RGE's as if all the scalar-bosons except h decouple at the scale M. 
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(c) If (71 (Q) < and (72 (Q) < 0, then 

+ _i_ y (MM _ (MM ."^j;^ _ > . (3.) 

If (7i(Q) < and X{Q) < 0, then 

^2 (Q) + ^ I ^{Xl{Q)X2{Q)-t{Q)) i^hlMpM-af^Q)^ _ (Q) ^(g) I > . (33) 

If (T2(Q) < and A(g) < 0, then 

^1 (0) + ^ I ^(Ai(Q)A2(Q)-A'(Q)) (^M^^^_^2(g)^ _ (Q) ^(g) | > q . (34) 

[ When (7i((5), (Q) and A((5) are all negative, the above three conditions are 
equivalent. ] 

In addition to the above conditions, we also demand local stability of the potential at the 
electroweak scale, namely, we calculate the mass spectrum of all scalar fields at the extremum 
of the potential and demand that all eigenvalues of the squared scalar mass are positive. We 
scan the remaining seven- dimensional space of Aj and cij and examine whether a given mass 
of the lightest CP-even Higgs boson is allowed under the above conditions. In this way we 
obtain the allowed range of as a function of tan P and M, for each value of the cut-off 
scale A. 

First, we discuss our result in the decoupling case, in which the soft-breaking scale M 
is much larger than the electroweak scale ~ v, and the masses of all the Higgs bosons but 
h (and 52) are at the order of M ^. In Fig. 2, the allowed range of ruh is shown as a 
function of tan/? for M — 1000 GeV. (We take the pole mass of top quark rrit — 175 GeV, 



^In the decoupling regime (M 00, which leads to a — > /? — ^ and X ^ 0); the masses of h and 
^2 are dominated by the (11) component of the mass matrix in Eq. (7) and the (22) component of 
that in Eq. (9), respectively. The mass of h is determined by the self-coupling constants A,, while 
that of S2 depends not only on the self-couplings constants ai but also on the free mass parameter 
mo. As noticed in the footnote 7, from ruh to M, the SM RGE are used in our analysis, even if 
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FIG. 2. The allowed mass range of the lightest CP-even Higgs boson for M = 1000 GeV. A is 
the cut-off scale. 

Oisi'iTT'z) — 0.118 for numerical calculation.) The allowed ranges are shown as contours for 
six different values of A, i.e. A = 10^^, 10l^l0^^10^^10^ and 10^ GeV. For most values 
of tan P, except for small tan /3 region, the upper bound of is about 175 GeV and the 
lower bound is between 110 GeV and 120 GeV for the cut-off scale A to be near the Planck 
scale. The numerical values in this figure are very close to those in the corresponding figure 
for the THDM discussed in Ref. [13]. Compared to the corresponding lower mass bound in 
the SM, which is 145 GeV when using the one-loop RGE's, the lower mass bound in this 
model is reduced by about 30 GeV to 40 GeV. The reason is similar to the THDM case: the 
lightest CP-even Higgs boson mass is essentially determined by the value of A2 for tan /3 to 



the mass of S2 is smaller than M. The effect of S2 on the mass bound of h is expected to be small, 
because at the one-loop level the primary effect is through the running of gi, whose contribution 
to the right-handed side of the RGE for the Higgs-self coupling constant is small. 
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FIG. 3. The allowed mass range of the lightest CP-even Higgs boson for M = 100 GeV. 

be larger than about 2 , where A2 plays the role of the self-coupling constant of the Higgs 
potential in the SM^. On the right-hand side of the RGE for A2, cf. Eq. (A5), there are 
additional positive-definite terms (A3 -|- (A3 -|- A4)^ -|- Ag -|- (7|) as compared to the RGE 
for the Higgs self-coupling constant in the SM. These additional terms can improve vacuum 
stability, and allow lower values of rrih. Therefore, one of the features of the model is to have 
a different mass range for the hghtest CP-even Higgs boson as compared to the SM Higgs 
boson, for a given cut-off scale. 

Next, we show our result for M to be around v. In Fig. 3, we present the ruh bound for 
M — 100 GeV. In this case, the allowed range of is reduced as compared to that in the 
decoupling case, and lies around mh ~ M for large tan /3. Notice that we have not included 



^However, tan /? cannot be too large to ignore the contribution of the bottom quark in the case 
with the type-11 Yukawa interaction. 
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phenomenological constraints from the b ^ s^, p parameter and the direct Higgs boson 
search experiment at LEP. As mentioned before, the mass bounds obtained from the RGE 
analysis are the same for the type-I and type-ll models without these phenomenological 
constraints. However, it was shown in Ref. [13] that the 6 — > 57 data can put a strong 
constraint on the allowed range of the Higgs boson mass for M < 200-400 GeV in the 
type-ll THDM, whereas there is no appreciable effect in the type-1 model. This is because a 
small M implies a light charged Higgs boson in the THDM which can induce a large decay 
branching ratio for 6 ^ S7 in the type-II model [19]^°. We expect a similar constraint from 
the 6 — > 57 data on the type-ll Zee-model, when M is small. 

In Fig. 4, we show the upper and lower bounds of as a function of M for various 
values of A. For given M, we scan the range of tan/3 for 1 < tan/3 < 16a/2 (~ 22.6). We 
find that the obtained m/j bounds are almost the same as those for the THDM. The primary 
reason for this is that the new coupling constants cti, (72 and (73 do not appear directly in 
the mass formula for m^, and therefore, do not induce large effects on the bounds of nih- 

We also investigate the allowed range of coupling constants ai, 02 and (J3. For this 
purpose, we fix a\ (or (72, (73) as well as tan/3 and M to evaluate the upper and the lower 
bounds of m/i for each A value. In this way, we determine the possible range of a\ (or (72, 
(73) under the condition that the theory does not break down below the cut-off scale A. In 
Fig. 5, we present the allowed range of Ox and for different choice of A in the case of 
M = 1000 GeV and tan /3 = a/2 or 16a/2. A similar fi gure is shown for the possible range of 
(T2 in Fig- 6. Wc see that the maximal value of o\ and 02 is around 0.7 for vrih = 110 — 170 
GeV if we take the cut-off scale to be 10^^ GeV. For smaller value of A the allowed ranges 
of Oi becomes larger. For example, a\ can exceed 1 for A = 10^^ GeV. We have calculated 
for other value of tan/3 and checked that these figures does not change greatly between 



In addition, it has been known that the i?b data also give strong constraints on the charged 
Higgs bosons in the type-II THDM [20]. 
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400 600 

M (GeV) 

FIG. 4. The allowed ranges of the lightest CP-even Higgs boson mass as a function of M for 
various A values. 

tan P — lA and 16-\/2. We also present the allowed range in the ai and (J2 plane for a fixed 
value of nih in Figs. 7 and 8 for mh = 125 GeV and rrih = 140 GeV, respectively. For either 
value of nifi with tan /3 = 16-\/2, both ai and (72 can be as large as 0.5 (2) for A = 10^^ (lO'^) 
GeV. The allowed range of (T3 and ruh for various values of A is given in Fig. 9. It is shown 
that, (73 has to be larger than zero, due to the vacuum stability condition. The maximal 
value of (73 is about 1 (3) for A = 10^^ (10^ ) GeV and M = 1000 GeV. The impact of these 
new coupling constants on the collider phenomenology is discussed in the next section. 
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FIG. 5. The allowed range of ai and ruh for various A values. 
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FIG. 6. The allowed range of (72 and ruh for various A values. 
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FIG. 7. The allowed range of ai and cr2 for rrih = 125 GeV. 
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FIG. 8. The allowed range of ai and cr2 for rrih = 140 GeV. 
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FIG. 9. The allowed range of as and ruh for various A values. 
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IV. TWO-PHOTON DECAY WIDTH OF THE NEUTRAL HIGGS BOSON 



In this section, we study the phenomenological consequences of the Higgs boson mass and 
the Higgs-boson-couphng constants derived in the previous section. The important feature 
of the Higgs sector of the Zee-model is that there are an additional weak doublet and a 
singlet charged Higgs boson. The physical states of the Higgs particles are two CP-even 
Higgs bosons, one CP-odd Higgs boson and two pairs of charged Higgs bosons. Therefore, 
the Higgs phenomenology is quite close to the ordinary two-Higgs-doublet model. One 
unique difference is the existence of the additional weak-singlet charged Higgs boson. The 
effect of this extra charged Higgs boson is especially important when M is much larger 
than the Z boson mass, i.e. in the decoupling regime. In such a case, the heavier CP-even 
Higgs boson, the CP-odd Higgs boson as well as one of the charged Higgs bosons have 
masses approximately equal to M, and these heavy states are decoupled from low energy 
observables. (Note that the condition on the applicability of the perturbation theory forbids 
too large self-couplings among the Higgs bosons. Hence, in the hmit of large M, the heavy 
Higgs bosons decouple from the low energy effective theory. ) The remaining light states 
are the lighter CP-even Higgs boson h and the lighter charged Higgs boson ^'2 which mainly 
comes from the weak-singlet. In the previous section, we show that, even in the decoupling 
case, there can be large difference in the allowed range of ruh between the Zee-model and the 
SM. Similarly, we expect that, even in the decoupling case, the presence of the additional 
weak-singlet charged Higgs boson can give rise to interesting Higgs phenomenology. 

Since the hghter charged Higgs boson 5*2 can couple to Higgs bosons and leptons, it can 
affect the decay and the production of the neutral Higgs bosons at colliders through radiative 
corrections. In the following, we consider the decay width of /i ^ 77 as an example. For 
a SM Higgs boson, the partial decay width (or branching ratio) of ^ 77 is small: ~ 9.2 
keV (or 2.2 x IQ-^) for rrih = 125 GeV, and - 15.4 keV (or 1.9 x IQ-^) for ruh = 140 GeV, 
with a 175 GeV top quark. Nevertheless, it is an important discovery mode of the Higgs 
boson at the LHC experiments for less than twice of the 1^-boson mass. Needless to 
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say that a change in the branching ratio of /i — > 77 would lead to a different production 
rate of pp hX — > 77X. At future e+e" LC's, the branching ratio of ^ 77 can be 
determined via the reaction e+e" —>■ qq^yj and e+e" —>■ 1/1777 with a 16-22% accuracy [16]. 
At the photon-photon coUision option of the future LC's, the partial decay width of /i — > 77 
can be precisely tested within a 2 % accuracy [17] by measuring the inclusive production 
rate of the Higgs boson h. Clearly, a change in the partial decay width of /i — > 77 will 
lead to a different production rate for h. In the Zee-model, such a change is expected after 
taking into account the loop contribution of the extra charged Higgs boson. We find that 
the deviation from the SM prediction can be sizable, and therefore testable at the LHC and 
future LC's. 

The partial decay width of ^ 77 is calculated at the one-loop order. Similar to our 
previous discussion, we limit ourselves to the parameter space in which 1 < tan/3 < 16\/2, 
and keep only the top quark contribution from the fermionic loop diagrams. Including the 
loop contributions from the W boson and the charged Higgs bosons and 5*2 together with 
the top quark loop contribution, we obtain [21] 



T{h 77) 



2567r2 sin Oy^m^ 



2 

--S-,,S2,t,W 



(35) 



with 



4 /cosa\ 
3 \smp J 

Iw = sin(/3 - a)Fi{ri) , 



where = and rrii is the mass of the internal lines in the loop diagram. Rs^ and Rs^ 
are given by 

^2 1 
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Rs, 



2 ml 



cos^ X |— Ai sin a sin^ P cos P + X2 cos a sin P cos^ P 



+X3 ^cosQ;sin^/3 — sinacos^/?^ — -(A4 -|- A5) cos(q; -|- P) sin2/3| 
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+ sin^ X sin a cos (5 -\- 02 cos a sin /3} + -\/2 sin^cosx— sin(Q; — j3) 



(36) 
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1 



sin^ X I — Ai sin a sin^ (3 cos /3 + A2 cos a sin (5 cos^ /3 



+A3 {cos a sin^ — sin o; cos^ ~ 9 (''^4 + ''^s) cos(q; + /?) sin 2/?| 



+ cos^ X {— Ci sin a; cos (3 -\- (72 cos a sin j3} — \/2 sinx cosx— sin(Q; — j3) 



(37) 



and 



Fo(r) = r (1 - r/(r)) , 
Fi/2(r) = -2r(l + (l-r)/(r)) , 
Fi(r) = 2 + 3r + 3r (2 - r) /(r) , 



(38) 
(39) 
(40) 



with 



f{r) 



sin ^ (^^Jl/rj 

lni+^ 



if r > 1 
if r < 1 



(41) 



l-Vl-r 

In the decouphng case of the model, namely ^ AjV^, the above formulae are greatly 
simplified. This limit corresponds to a ^ /3 — | and x — 0, so that the hght charged Higgs 
boson S2 is identical to the weak-singlet Higgs boson u^. Thus, we have 

1 



R 



'S2 



(cTi cos^ P + (72 sin^ , 



(42) 



and both the top-quark and the W boson loop contributions reduce to their SM values. We 
like to stress that the weak-singlet Higgs boson does not directly couple to the quark fields 
in the hmit of x — 0. Therefore, it does not affect the decay rate of 6 — > 57 at the one-loop 
order. Similarly, being a weak singlet, it also gives no contribution to the p parameter. 
Hence the low-energy constraint from either the h ^ decay or the p parameter on the 
Zee-model in the limit of x — is similar to effects of that on the THDM. Let us examine 
at the one-loop effect of the weak-singlet charged Higgs boson on the decay width of /i — 77 
in the decoupling limit. Let us recall that in Fig. 8, the size of the new couplings (7i and 
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FIG. 10. (a) The ratio r as a function of the charged Higgs boson mass ms^ for various values 
of the couphng constants a\ = 02 = o with = 125 GeV. The two smaller cr's are consistent with 
the cut-off scales A = 10^^ GeV and A = 10^^ GeV, respectively. The two larger cr's are allowed 
for A = 10"^ GeV. (b) A similar plot with m/j = 140 GeV. 

(72 can be as large as 2 simultaneously, if the cut-off scale is at the order of 10^ GeV. For 
the Zee- model to be a valid low energy effective theory up to 10^^ GeV, a\ and 02 cannot 
be much larger than 0.6. To illustrate the imphcations of this result, we show in Figs. 10 
(a) and 10 (b) the ratio (r) of the /i ^ 77 width predicted in the Zee-model to that in the 
SM, r = rzee(^ Jj)/^SM{h 77), as a function of the coupling constant (72 and the 
charged Higgs boson mass . Here, for simplicity, we have set ai — so that the tan /3 
dependence drops in the decoupling case, cf. Eq. (42). For illustrations, we consider two 
cases for the mass of the lighter CP-even Higgs boson: — 125 GeV and rrih = 140 GeV. 



As shown in the figures, the ratio r can be around 0.8 for iTi = (72 = cr ~ 0.5 and 



100 



GeV. This reduction is due to the cancellation between the contribution from the -So-boson 
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FIG. 11. (a) The ratio r as a function of the charged Higgs boson mass for negative values 
of the coupUng constants (T2 with ruh = 125 GeV, ai = and tan/3 = 16\/2. The value a2 = —0.2, 
—0.5 or —0.8 is consistent with the cut-off scale A = 10^^, 10^ or 10^ GeV, respectively, (b) A 
similar plot with m^j = 140 GeV, ci = and tan/3 = 16\/2- The value (T2 = —0.25, —0.6 or —1 is 
consistent with the cut-off scale A = 10^^, 10^ or 10^ GeV, respectively. 

loop and the W^-boson loop contributions. To have a similar reduction rate in Tzeeih 77) 
for a heavier S2, the coupling constant (72 (and cti) has to be larger. Next, as shown in 
Figs. 7 and 8, ai and (T2 do not have to take the same values in general, and they can be less 
than zero. In the case where both ai and (T2 are negative, the contribution of the S'2-loop 
diagram and that of the VF-loop diagram have the same sign, so that r can be larger than 
1. Such an example is shown in Fig. 11 (a), where the ratio r for ruh — 125 GeV is shown as 
a function of rrig^ at various negative (T2 values with ai = and tan P = 16\/2. We consider 
the case with a2 — —0.2, —0.5 or —0.8, which is consistent with the cut-off scale A = 10^^, 
10^ or 10^ GeV, respectively. In the case of A = 10^^ GeV (10^ GeV), the deviation from the 
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FIG. 12. (a) The partial decay width T{h S^S^) for ms^ = 80,100,150,200 GeV with 
0-1 = 0-2 = 1 for the allowed range of from 100 GeV to 500 GeV. (b) The ratio of T{h S^S^) 
with the total decay width of the SM Higgs boson for each values of m^^- 

SM prediction can be about +6% (+30%) for rrig^ = 100 GeV. In Fig. 11 (b), the similar 
plot of the ratio r is shown for nih — 140 GeV with ai — and tan /3 = IdV^. Each case 
with (72 — —0.25, —0.6 or —1 is consistent with A = 10^^, 10^ or 10^ GeV, respectively. The 
correction is larger in the case with nih = 140 GeV than in the case with mh = 125 GeV 
for a given A. The deviation from the SM prediction can amount to about +8% (+40%) for 
A = 10^^ GeV (10^ GeV) when rrig^ — 100 GeV. Larger positive corrections are obtained 
for smaller rrig^ values. Such a deviation from the SM prediction can be tested at the LHC, 
the e+e~ LC and the 77 option of LC. 
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Before concluding this section, we remark that if nih is larger than 2ms2 such that the 
decay mode h —>■ S2S2 is open, the total decay width of h can be largely modified from the 
SM prediction for large ai^2- In terms of Rs2, the partial decay width of /i — > S2S2 is given 

by 



r(h ^ 8^82) = 



'"^^^ (43) 



2 ' 
mi 



IGnmh ^ 

where = {2m%^Rs2/v'^Y . In Fig. 12 (a), we show the partial decay width V{h 82 82) 
for mg^ = 80, 100, 150, 200 GeV with ai = era = 1, cf. Eq. (42), for the allowed range of mn 
from 100 GeV to 500 GeV. In Fig. 12 (b), the ratio of r{h 8^82) to the total width of the 
SM Higgs boson (r^°*^'(5'M)) is shown as a function of for each value of m^^. This is to 
illustrate the possible size of the difference between the total width of the lightest CP-even 
Higgs boson h in the Zee-model and that of the SM Higgs boson^^. Clearly, the impact of 
the 82 82 decay channel is especially large in the small mh region. We note that r^°*^'(5'M) 
can be determined to the accuracy of 10-20% at the LHC and the LC if mh < 2m^, and to 
that of a couple of per cents if > 2771^ [23]. (m^ is the mass of the Z-boson.) Hence, 
measuring the total width of the lightest neutral Higgs boson can provide a further test 
of the Zee- model for m^ > ^mg^. The change in the total width also modifies the decay 
branching ratio of /i ^ ZZ, hence yields a different rate of /i ^ ZZ fi^ ^~ ^'^ ii~ for a 
given mh- (In the SM, the branching ratio of /i — > ZZ is about 1/3 for mh > 200 GeV.) 
Needless to say that for mh > ^mg^-, the production mode olh ^ 82 82 £'^i''^T is also 
useful to test the Zee-model. Further discussion on this possibility will be given in Sec. VI. 



^■"^In doing this analysis, we have in mind a low cut-ofF scale A = 10^ GeV, which allows a wide 
range of values for a's, mg^ and nih- 
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V. PHENOMENOLOGY OF CHARGED-HIGGS BOSONS 



In the Zee-model, two kinds of charged Higgs bosons appear. If there is no mixing 
between them (x = 0), the mass eigenstates Sf" and correspond to the THDM-hke 
charged Higgs field and the singlet Higgs field o;^, respectively. The case with x = occurs 
in the hmit of » i"^,/^^ and ml; i.e. in the decoupling hmit. The detection of St can 
be a clear indication of the Zee-model. As to be shown later, its phenomenology is found to 
be drastically different from that of the THDM-like charged Higgs bosons Sf [24] . Here, we 
discuss how the effects of this extra charged boson can be explored experimentally. We first 
consider the case with x = 0, and then extend the discussion to the case with a non-zero x- 

The S2 boson decays into a lepton pair e^I^^. with the coupling constant fij. The partial 
decay rate, Vfj — r(5'^ — > e'^T'l.), is calculated as 

and the total decay width of S2 is given by 

r§tai=Ert. (45) 

By taking into account the hierarchy pattern of /y, cf. Eqs. (23) and (24), and by assuming 

nig^ = 100 GeV and I/12I = 3 x 10"'', the total decay width and the life time (r) is estimated 



to be 



12 



r§tai~rf| + rf^ ~i.6kev, (46) 

r ~ l/r^^^ai ~ lO-^'sec. (47) 
This implies that 5*2 decays after travehng a distance of ~ 10~^° m, which is significantly 



^^The size of the decay width depends on the value of /12. If we take ms^ > 500 GeV or ^ < 100 
GeV, /12 can become one order of magnitude larger than 3 x 10~^, while still being consistent with 
the phenomenological bounds discussed in Sec. II. 
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shorter than the typical detector scale. Therefore, S2 decays promptly after its production, 
and can be detected at collider experiments. 

The main production channel at the LEP-II experiment may be the pair production 
process e+e~ — > 3282, similar to the production of the THDM-like charged Higgs boson 
Si - The matrix-element squares for the S^'Sf production (i = 1, 2) are given by 

|A^(eZ(^)ei.) - Srsn\' = " ^(ll - .^^.Q^J^^^^^^^P^/?! sin^ ©, (48) 

where Qe = —1 and II = —\ (0) for the incoming electron (e^); Qs = ~1 ^^'^ ^s~~\ 



(0) for i = 1 (2); (3s. — -y^l — 4m|./s, sw — sin^vF, cw — cosOw, and © is the scattering 
angle of 5"" in the e^e~ center-of-mass (CM) frame whose energy is ^/s. For the other 
electron-positron helicity configuration (e^e^ and e^^e^), the cross sections are zero. Thus 
the total cross section for the S2S2 pair production is given by 



^ ^ 967r ^2 



1 



■ (49) 



Hence, the production rates of Si and S2 are different. We note that the ratio of cross 
sections for SiS^ and S2S2 production, a{e'^e~ S'2'5'^)/cr(e"'"e~ —>■ Si Si), is 0.8 at 
a/s = 210 GeV assuming that the masses of Sf and 5"^ are the same. This ratio is inde- 
pendent of the masses of Si and 5*2 for a fixed CM energy. (Only the difference between 
Si Si Z and S2 S2 Z coupling constants determines this ratio. ) 

The lower mass bound of the THDM-like charged boson Sf can be obtained by studying 
its TV and cs decay modes, completely in the same way as the charged Higgs boson search in 
the minimal supersymmetric standard model (MSSM) [22] . Similar experimental constraints 
may be obtained for the extra charged bosons Sf. The situation, however, turns out to be 
fairly different from the Sf case. First of all, decays of 5"^ are all leptonic. Secondly, the 
branching ratios of various 5"^ decay modes are estimated as 

5(^2" ^ e^^T) ~ 0.5, (50) 

B{S2 ^ H~$t) ^0.5, (51) 

-f j ~ 10-^ (52) 
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where we have used the relations given in Eqs. (23) and (24). Clearly, the branching ratio 
into the t'I^t mode is very small, so that it is not useful for detecting Sf at all. This is 
different from the case of detecting the ordinary THDM-like charged Higgs boson, which 
preferentially decays into heavy fermion pairs (e.g. tu and cs). Instead of studying the 
T^v'^ mode, the v'^ and //^ modes can provide a strong constraint on the mass of 
Sf. In fact, the branching ratio of S2 e~^T or h'^t is almost 100 %, so that we have 
a{e^e~ — > S2S2 — > t^l'^^r) ~ cr(e'*"e~ — > S2S2), where and represent e~ or //~ (not 
t"). Let us compare this with the cross section a{e^e'' W'^W' £'^i'~l^T) = cT(e^e~ — > 
W+W-)-B{W- t-^T?, where BiW' t'^T) = B{y^- e-^T)+B{W- i^-^t) - 
21%. As seen in Fig. 13, the cross section (7(e+e~ — > S2S2 — > 1'^1'~^t) is comparable with 
a{e+e- W+W- £+f-^T )• Therefore, by examining the LEP-II data for £+£'-^t 
(where i^i'^ = e+e^, e^/x^ or jJb^jjT ., in contrast to t+t~ for the Sf case), the experimental 
lower bound on the mass of S2 can be determined. Such a bound can be induced from 
the smuon search results at the LEP experiments [25,26] in the case that neutralinos are 
assumed to be massless. From the fi'^fi'I^T data accumulated up to \/s = 202 GeV [26], 
we find that the lower mass bound of is hkely to be 80-85 GeV for the x = cases. [ 
We note that the right-handed smuon 01%) in the MSSM carries the same SU{2) x U{1) 
quantum number as the weak-singlet charged Higgs boson {S2 for ~ 0). ] 

We next comment on 5'^-production processes at hadron coUiders and future LC's. At 
hadron colliders, the dominant production mode is the pair production through the Drell- 
Yan-type process. The cross sections for pp S2S2 at the Tcvatron Run-II energy [^/s = 2 
TeV) and pp — > S2S2 at the LHC energy {^/s — 14 TeV) are shown as a function of 
in Fig. 14 for x = 0. At future LC's, the S2 boson may be discovered through the above- 
discussed pair-production process from the electron-positron annihilation if y/s/2 > ms2- In 
Fig. 15, we show the total cross section of e"'"e~ S2 S2 for x = as a function of ms^ for 

= 300, 500, and 1000 GeV. 

Finally, we like to discuss the case with a non-zero X; in which S2 is a mixture of the 
singlet charged Higgs boson state (a;~) and the doublet charged Higgs boson state {H~). Let 
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FIG. 13. The cross section of the leptonic decay process e^e - 
i and ^' = e or ^) at ^ = 190, 200, 210 GeV. The process e+e" 
-y/s = 210GeV is shown for comparison. 

us see how the above discussion is changed in this case. The doublet charged Higgs bosons 
with the mass of 100 GeV mainly decay into the t~v and cs channels. Thus, the branching 
ratio of the decay process — > t'^T , where £~ represents e~ and is expressed in a 
non-zero x case as 









X=0 + C0S2 X 





(53) 



where r^otailx=o (^ — 1)2) is the total width of at x = with the same mass as the 
decaying S2 on the left-hand side of the above equation. The formula of rj^tailx=o is given 
in Eq. (45), while r^\3^i|^=o, which is the same as the total decay width of the charged Higgs 



boson in the THDM is given by 



totallx=0 - ^//')' 
If 



(54) 



where //' are fermion pairs which are kinematically allowed. In the type-11 Yukawa cou- 
plings, we have 
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FIG. 14. The total cross sections of pp — S2S2 at s/s = 2 TeV (solid curve) and pp — S2S2 



at ^/s = 14 TeV (dotted curve) as a function of nig^ 
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FIG. 15. The total cross section of e^e~ — S2S2 as a function of m^^ at ^/s = 300, 500 and 
1000 GeV. 
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FIG. 16. The decay branching ratio of S2 — t'^r (where ^~ = e~ or fi~) as a function of the 
mixing angle x for iT^Sn — GeV, tan/3 = 1 and various values of the coupling constant /12. 
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In the THDM, the total decay width of the charged Higgs boson {H~) for tuh- = 100 GeV 
is about 470 keV. Hence, if the mixing angle x is not so small, the decay pattern of 5"^ is 
dominated by that of the THDM charged Higgs boson H~. In Fig. 16, we plot the branching 
ratio B{S2 —>■ i~^T) as a function of sinx at m^^ = 100 GeV for several values of /12. We 
only show the case with tan (3 — 1, where the result is independent of the type of the Yukawa 
interaction. The couphng constant /12 is taken to be 3, 9, 18 and 36 (xlO~^), which satisfy 
the phenomenological constraints given in Sec. II. As expected, the branching ratio decreases 
as X increases. When /12 = 36 x 10~^, -6(5'^ i'^r) is smaller than 10% for sinx > 0.89. 
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For the smaller /12 values, the branching ratio reduces more quickly. The branching ratio is 
not sensitive to mg^ unless the mass exceeds the threshold of the decay into a th or h^W" 
pair. Above the threshold of the tb pair production, the decay rate of S2 tb is large due 
to the large mass of the top quarks, so that -6(5'^ — > £~^t) is substantial only for very 
small values of x- Finally, while the decay branching ratio can change drastically depending 
on the mixing angle x, the production cross section for e+e" — > S2S2 remains unchanged. 
In conclusion, the process e+e~ — > S2S2 — > £'^£'~^t can also be useful for testing the 
Zee- model in the non-zero x case, provided sin x is not too large. 

VI. DISCUSSION AND CONCLUSION 

In this paper, the Higgs sector of the Zee-model has been investigated, in which neutrino 
masses are generated radiatively. This model contains extra weak-doublet Higgs field and 
singlet charged Higgs field. 

We have studied indirect effects of these extra Higgs bosons on the theoretical mass 
bounds of the lightest CP-even Higgs boson, which are obtained from the requirement that 
the running coupling constants neither blow up to a very large value nor fall down to a 
negative value, up to a high-energy cut-off scale A. For A = 10^^ GeV, the upper bound 
of ruh is found to be about 175 GeV, which is almost the same value as the SM prediction. 
In the decoupling regime (M ^ ^2), the lower bound is found to be about 100 GeV for 
A = 10^^ GeV, which is much smaller than the lower bound in the SM, and is almost the 
same as that in the THDM. For smaller A values, the bounds are more relaxed, similar 
to that of the SM. We have also investigated the allowed range of the coupling constants 
relevant to the weak-singlet Higgs field. 

The most striking feature of the Zee-model Higgs sector is the existence of the weak- 
singlet charged Higgs boson. We have examined the possible impact of the singlet charged- 
Higgs boson on the neutral Higgs boson search through radiative corrections. We found 
that its one-loop contributions to the /i — > 77 width can be sizable. In the allowed range 
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of the coupling constants the deviation from the SM prediction for this decay width can be 
about —20% or near +10% for nig^ = 100 GeV and A = 10^^ GeV, depending on the sign 
of the couphng constants CTj. The magnitude of the deviation is larger for lower A values 
or for smaller ms^ values. For example, a positive deviation over 30-40% is possible for 
= 125-140 GeV, m^^ = 100 GeV, and A = 10^ GeV. 

In the decoupling limit (i.e. when ^ v"^, where a — > /3 — 7r/2 and x ~^ 0); we 
expect that the production cross sections for gg — > h, e^e~ — > vVh and e^e' — > Z^h in 
the Zee-model are the same as those in the SM. However, a sizable change in the decay 
branching ratio of /i — >^ 77 can alter the production rate of pp — > hX — > 77X at the 
LHC, where this production rate can be determined with a relative error of 10-15% [15]. 
Also, such a deviation in the branching ratio of /i ^ 77 directly affects the cross section 
of e+e" — > vVh (andZ^/i) — > 1/1777, which can be measured with an accuracy of 16-22% 
at the future e+e~ LC (with ^/s — 500 GeV and the integrated luminosity of 1 ab"^) [16]. 
Therefore, the Zee-model with low cut-off scales can be tested through the ^ 77 process at 
the LHC and the e+e^ LC's. At the future photon colliders, the enhancement (or reduction) 
of the /i — > 77 partial decay rate will manifest itself in the different production rate of h 
from the SM prediction. A few per cent of the deviation in T{h — > 77) • B(h ^ bb) can be 
detected at a photon collider [17], so that the effects of the singlet charged Higgs boson can 
be tested even if the cutoff scale A is at the Planck scale. 

The collider phenomenology of the singlet charged Higgs boson has turned out to be 
completely different from that of the THDM-likc charged Higgs boson. The singlet charged 
Higgs boson mainly decays into £^^t (with = or while the decay mode t^^t 
is almost negligible due to the relation I/12I » I/13I » I/23I. This hierarchy among the 
coupling constants fij results from demanding bi-maximal mixings in the neutrino mass 
matrix generated in the Zee-model to be consistent with the neutrino oscillation data. On 
the other hand, the THDM-hke charged Higgs boson decays mainly into either the tv mode 
or the cs mode, through the usual Yukawa-interactions. Hence, to probe this singlet charged 
Higgs boson using the LEP-II data, experimentalists should examine their data sample with 
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e^e~^T , e'*'/i~^r , A*"'"e~^T or ^^^~^t , while the experimental lower mass bound of the 
THDM-like charged Higgs boson is obtained from examining the tt^t , t^t jj and jjjj 
events. Using the published LEP-II constraints on the MSSM smuon production (assuming 
the lightest neutralinos to be massless), we estimate the current lower mass bound of this 
singlet charged Higgs boson to be about 80-85 GeV. The Tevatron Run-II, LHC and future 
LC's can further test this model. 

Finally, we comment on a case in which the singlet charged Higgs boson (S^ for x = 0) 
is the hghtest of all the Higgs bosons. For mh/2 > > rnz, the Higgs sector of the 
Zee-model can be further tested by measuring the production rate of pp (or pp) —>■ hX 
S}S2X ^+f-^r X. The branching ratio for h S}S2 ^+^'~^t can be large. For 
instance, for rrif^ = 210 GeV and rrig^ — 100 GeV, this branching ratio is about 12% for 
each = e+e~, ii~ , ii^e~ or ii^ ii~ . The branching ratio decreases for larger masses 

of h. Moreover, the total decay width of h can be largely modified when the decay channel 
h — > S2S2 is open. In this case, the decay branching ratios of /i — > ZZ are also 

different from the SM predictions. 

In conclusion, the distinguishable features of the Zee-model from the SM and the THDM 
can be tested by the data from LEP-II, the Tevatron Run-II and future experiments at LHC 
and LC's. 

ACKNOWLEDGMENTS 

We are grateful to the warm hospitality of the Center for Theoretical Sciences in Taiwan 
where part of this work was completed. CPY would like to thank H.-J. He, J. Ng and 
W. Repko for stimulating discussions. SK was supported, in part, by the Alexander von 
Humboldt Foundation. GLL and JJT were supported, in part, by the National Science 
Council of R.O.C. under the Grant No NSC-89-2112-M-009-035; YO was supported by the 
Grant-in-Aid of the Ministry of Education, Science, Sports and Culture, Government of 
Japan (No. 09640381), Priority area "Supersymmetry and Unified Theory of Elementary 



39 



Particles" (No. 707), and "Physics of CP Violation" (No. 09246105); CPY is supported by 
the National Science Foundation in the USA under the grant PHY-9802564. 



40 



APPENDIX: ONE-LOOP RGE'S FOR DIMENSIONLESS COUPLING 

CONSTANTS 

Here, we summarize the relevant RGE's to our study. For the gauge couphng constants, 
we have 

(Al) 
(A2) 
(A3) 

The RGE's for the Higgs-self-couphng constants of the doublets are calculated at the one- 
loop level as 
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1 22 
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/^T-^i = |l2A? + AXl + 4A3A4 + 2X1 + 2A^ + 2a? 
d^ lon^ ^ 



- {Sgl + 9gl) A, + [^gt + \glgl + lat) } , (A4) 
1^-^X2 = {12A2 + AXl + 4A3A4 + 2X1 + 2A^ + 2al + 12y^X2 - 

- {39! + 9,,^) X2 + [lat + \9l9l + Igl) ] , (A5) 

/i|^A3 = {2 (Al + A2) (3A3 + A4) + AXl + 2Al + 2A^ + 2(7i(72 + Gy^Xs 

- {Sgj + 9gl) A3 + [Igt - \glgl + \gi) ] , (A6) 

1^-^X4 = {2 (Al + A2) A4 + 4 (2A3 + A4) A4 + 8X1 + %t A4 

-{^9l + Q9l)>^4 + 3glgi}, (A7) 
/i^As = ^ {2A1 + 2A2 + 8A3 + I2A4 + - (3^2 ^ 9^2^ I (A8) 

and those with respect to the additional singlet charged Higgs are given by 
d 1 



{Aal + 2(71(73 + 6Ai(7i + (4A3 + 2A4) (72 + Sfijfijai 



- {y^I + l9i) <Ti + 3^^} , (A9) 
Ai^o-2 = {"^^2 + 2cr2a3 + 6A2cr2 + (4A3 + 2A4) cxi + 6y^a2 + 8fijfija2 
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y^7? + Igl) ^2 + 3,^} , (AlO) 
/^^fT3 = ^ {Sal + 8al + hal + 16/,,^(73 - 128 tr f - \2g\a^ + 24^^} . (All) 

Finally, the RGB's for the Yukawa-type couphng constants are obtained at one-loop level as 
'^Ifl^'^ ^ {~ (^^^ ^ ^^^) -^'^ ^ ^fkifkifij - 4/ifc/w//j| , (A13) 

where 

tl^ = XI fijfjkfklfli, 
i,j,k,l=l—3 

fijfij = X fijfij- 
i,j=l-3 
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